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ABSTRACT 

The theoretical basis for the link between the leaf exchange 
of carbonyl sulfide (COS), carbon dioxide (CO2) and water 
vapour (H 2 0) and the assumptions that need to be made in 
order to use COS as a tracer for canopy net photosynthesis, 
transpiration and stomatal conductance, are reviewed. The 
ratios of COS to CO2 and H 2 0 deposition velocities used to 
this end are shown to vary with the ratio of the internal to 
ambient CO> and H 2 0 mole fractions and the relative limi- 
tations by boundary layer, stomatal and internal conduc- 
tance for COS. It is suggested that these deposition velocity 
ratios exhibit considerable variability, a finding that chal- 
lenges current parameterizations, which treat these as 
vegetation-specific constants. COS is shown to represent a 
better tracer for CO2 than H 2 0. Using COS as a tracer for 
stomatal conductance is hampered by our present poor 
understanding of the leaf internal conductance to COS. 
Estimating canopy level C0 2 and H 2 0 fluxes requires dis- 
entangling leaf COS exchange from other ecosystem 
sources/sinks of COS. We conclude that future priorities for 
COS research should e to improve the quantitative under- 
standing of the variability in the ratios of COS to C0 2 and 
H 2 0 deposition velocities and the controlling factors, and to 
develop operational methods for disentangling ecosystem 
COS exchange into contributions by leaves and other 
sources/sinks. To this end, integrated studies, which concur- 
rently quantify the ecosystem-scale C0 2 , H 2 0 and COS 
exchange and the corresponding component fluxes, are 
urgently needed. 
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Abbreviations: a, coefficient in Eqn 4; CA, carbonic anhy- 
drase; CO2, carbon dioxide; COS, carbonyl sulfide; D x , dif- 
fusivity of C0 2 (D c ), H 2 0 (7> v ) or COS (7> s ) in air; e x , 
ambient (e„) or internal (ej) water vapour pressure (kPa); 
_F e s , ecosystem-scale flux of COS (pmol m~ 2 s _1 ); 7i x , leaf flux 
of CO2 (Ff, jUtnol nr 2 s" 1 ), H 2 0 (7i v , mmol nr 2 s" 1 ) or COS 
(71 s , pmol m~ 2 s -1 ); 7ni s , non-leaf above-ground flux of COS 
(pmol nr 2 s" 1 ); 7 S S , soil flux of COS (pmol m" 2 s" 1 ); G b s , 
combined aerodynamic and quasi-laminar boundary layer 
conductance for COS (mol m" 2 s" 1 ); g b x , leaf boundary layer 
conductance for C0 2 (g b c , mol nr 2 s -1 ), H 2 0 (g b v , mol 
m~ 2 s -1 ) or COS (g b s , mol nr 2 s _1 ); g c h C , chloroplast conduc- 
tance for CO2 (mol m" 2 s -1 );g gp x , gas phase conductance for 
C0 2 (ggp C , mol m~ 2 s -1 ) or COS (ggp S , mol m" 2 s" 1 ) in the 
intercellular space; Gi S , canopy-scale internal conductance 
for COS (mol m" 2 s -1 ); g x , leaf internal conductance to C0 2 
(gi C , mol nr 2 s -1 ) or COS (gi S , mol m~ 2 s _1 ); GPP, gross pho- 
tosynthesis (/itnol in" 2 s _1 ); G s s , canopy-scale stomatal 
conductance for COS (mol m" 2 s" 1 ); g s x , leaf stomatal con- 
ductance to C0 2 (g s c , mol m" 2 s" 1 ), H 2 0 (g s v , mol m" 2 s" 1 ) or 
COS (g s s , mol nr 2 s -1 ); g w x , combined cell wall, plasma mem- 
brane and cytosol conductance for CO2 (g w c , mol m~ 2 s _1 ) or 
COS (g w s , mol m" 2 s _1 ); H 2 0, water vapour; H 2 S, hydrogen 
sulfide; NEE, net ecosystem CO2 exchange (,umol m" 2 s" 1 ); 
Rhs, ratio of leaf boundary layer to stomatal conductance for 
COS; RECO, ecosystem respiration (jUmol m" 2 s" 1 ); RH, 
relative humidity (fraction); R si , ratio of leaf stomatal to 
internal conductance for COS; s, slope of water vapour 
saturation function (°C _1 ); 7 a , air temperature (°C); 71, leaf 
temperature (°C); VPD, vapour pressure deficit (kPa); f3 x , 
beta coefficient for CO2 (If) or H 2 0 (/3 V ); A7, difference 
between air and leaf temperature (°C); X x , ratio of deposi- 
tion velocities of COS to CO2 (A c ) or COS to H 2 0 (A v ); #, x , 
ambient CO2 (% a c , p.mo\ mol" 1 ), H 2 0 (^ a v , mmol 
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mol -1 ) or COS (j£a S , pmol mol -1 ) mole fraction; % c c , chloro- 
plast CO2 (jUmol mol -1 ) mole fraction; j£ cs x , chloroplast 
surface C0 2 (%cs C , jUmol mol -1 ) or COS Gfcs S , pmol mol -1 ) 
mole fraction; # x , CO2 Gfr c , //mol mol -1 ), H 2 0 (# v , mmol 
mol -1 ) or COS (# s , pmol mol -1 ) mole fraction in the sub- 
stomatal cavity; ^ w x , cell wall CO2 (%w C , jUmol mol" 1 ) or COS 
CtV s , pmol mol -1 ) mole fraction. 

INTRODUCTION 

Canopy net photosynthesis, transpiration and stomatal con- 
ductance are key conceptual terms in most contemporary 
models of ecosystem carbon and water cycling (Sitch et al. 
2008). While net photosynthesis, transpiration and stomatal 
conductance can be quantified accurately with enclosures at 
the leaf scale [e.g. von Caemmerer & Farquhar (1981), but 
see Rodeghiero, Niinemets & Cescatti (2007)], obtaining 
reliable estimates at the canopy scale is much more difficult. 
Scaling up net photosynthesis, transpiration and stomatal 
conductance measured in leaf enclosures to the canopy 
scale requires minimum knowledge on (1) the response of 
leaf gas exchange rates to environmental drivers; (2) how 
this response changes with depth in the canopy; (3) the 
vertical variation of environmental drivers within the plant 
canopy; and (4) the vertical distribution of the assimilating/ 
transpiring plant area (Kruijt, Ongeri & Jarvis 1997; Wohl- 
fahrt et al. 2010). Transpiration of individual trees can be 
directly quantified by sap flux methods; however, up-scaling 
methods are again needed to turn these measurements into 
canopy transpiration (Wilson et al. 2001). Enclosures that 
include entire ecosystem are prone to artefacts due to 
modifications of the environmental conditions (Dore et al. 
2003). In addition, ecosystem enclosures yield the net eco- 
system CO2 and H 2 0 exchange, that is, they are unable to 
partition between canopy net photosynthesis/transpiration 
and CO2/H2O fluxes from/to the soil and other ecosystem 
components. This drawback also holds for micrometeoro- 
logical techniques such as the eddy covariance method, 
which, however has the advantage of being unobtrusive and 
able to provide near-continuous long-term flux data (Bal- 
docchi, Hicks & Meyers 1988; Baldocchi 2003). In order to 
recover canopy net photosynthesis and transpiration from 
net ecosystem fluxes of CO2 and H 2 0, it is necessary to 
concurrently quantify/estimate the confounding CO2 and 
H 2 0 fluxes - a non-trivial task in particular for C0 2 , which 
has multiple sources in an ecosystem. Converting estimates 
of canopy transpiration to stomatal conductance in turn is 
fraught with problems due to difficulties with the correct 
specification of the vapour gradient between the transpiring 
surface and ambient air (Magnani et al. 1998). As a conse- 
quence, available estimates of canopy net photosynthesis, 
transpiration and stomatal conductance are inherently 
uncertain, which in turn translates into uncertain model 
parameterizations and predictions. 

Recently, several authors have advocated measurements 
of COS exchange to provide independent constraints on 
canopy net photosynthesis in particular (Sandoval-Soto 
etal. 2005; Montzka et al. 2007; Campbell etal. 2008; 



Brugnoli & Calfapietra 2010; Seibt et al. 2010; Stimler et al. 
2010a, 2011), as well as on canopy transpiration and sto- 
matal conductance (Seibt et al. 2010). The rationale for 
these proposals derives from both leaf and (very few) eco- 
system flux measurements, which show a high degree of 
co-variation between the net exchange rates of CO2, H 2 0 
and COS (Xu, Bingemer & Schmidt 2002; Sandoval-Soto 
et al. 2005; Stimler et al. 2010a). 

Given the promising possibility of quantifying 
ecosystem-scale COS exchange by using the eddy covari- 
ance method in combination with new analytical techniques 
(Graus, Muller & Hansel 2010; Stimler, Nelson & Yakir 
2010b), in order to better constrain canopy photosynthesis, 
transpiration and stomatal conductance, the objective of the 
present paper is: (1) to review the mechanistic link between 
leaf- and ecosystem-scale CO2, H 2 0 and COS fluxes; (2) to 
critically evaluate the assumptions required for estimating 
canopy net photosynthesis, transpiration and stomatal con- 
ductance from COS exchange measurements; and finally 
(3) to indicate areas of future research. 

THE LINK BETWEEN LEAF C0 2 , H 2 0 AND 
COS EXCHANGE 

We begin our assessment with a review of the equations 
describing the diffusive flux of CO2, H 2 0 and COS in/out of 
leaves, as sketched in Fig. 1. Leaf net photosynthesis (Ff) is 
given as 




Figure 1. Schematic illustration of the diffusion pathways of 
CO2, H2O and COS into/out of a leaf. Blue panels represent mole 
fractions; red panels are conductances using abbreviations 
introduced in the text. Additional abbreviations include: g gp c and 
g gp s referring to gas phase conductances for CO2 and COS in the 
intercellular space; g w c and g w s (including cell wall, plasma 
membrane and cytosol) and % w c an d £„ s referring to cell wall 
conductances and mole fractions of CO2 and COS, respectively; 
£ cs c and £ cs s referring to chloroplast surface mole fractions of 
CO2 and COS; g C h C referring to the chloroplast conductance to 
CO2 and Xc C referring to the chloroplast CO2 mole fraction. 
Following Stimler et al. (2010a), who concluded CA to be 
effectively located at the chloroplast surface, the end point for 
diffusion of COS is placed at the inside of the chloroplast 
surface. 
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and the leaf net flux of COS (Ff) as 
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Here, % refers to mole fractions (subscripts: a . . . ambient 
air, i . . . intercellular space, c . . . chloroplast) and g to con- 
ductances (subscripts: b . . . boundary layer, s . . . stomata, 
i . . . internal). Note that we employ a sign convention by 
which fluxes directed into the leaf have a negative sign. The 
boundary layer and stomatal conductances may be inter- 
converted between COS, CO2 and H 2 0 based on their dif- 
fusivity (D) ratios, that is, 



(4) 



where the exponent 'a' takes a value of 1 for g s (i.e. molecu- 
lar diffusion) and 0.67 for g h (i.e. forced convection which 
typically applies in well-ventilated leaf enclosures) (Camp- 
bell & Norman 1998). g s s and gb S may thus be converted to 
their CO2 (H 2 0) counterparts by multiplication with ca. 
1.21 (2.00) and 1.14 (1.59), respectively (Seibt etal. 2010; 
Stimler et al. 2010a). 

While Eqns 1 and 2 are well established (von Caemmerer 
& Farquhar 1981), Eqn 3 merits further explanation: the 
endpoint of the diffusion gradient for COS is the location of 
the enzyme CA, which appears to be available throughout 
the plasma membrane, cytosol, chloroplast envelope and 
stroma (Evans et al. 2009). Stimler etal. (2010a), following 
Gillon & Yakir (2000), concluded CA to be effectively 
located at the chloroplast surface, an assumption we graphi- 
cally followed in Fig. 1. CA, which has an extremely high 
sensitivity to COS [larger by a factor of 1000 compared with 
C0 2 ; Protoschill-Krebs, Wilhelm & Kesselmeier (1996)], is 
responsible for the hydration of COS, as nicely demon- 
strated by Stimler et al. (2011) with CA-deficient antisense 
lines of C3 and C4 plants. The hydration of COS is essentially 
a one-way reaction whereby one molecule H 2 S and C0 2 are 
generated for each hydrated molecule COS (Protoschill- 
Krebs et al. 1996). As emissions of COS out of leaves have 
not been reported even at very low ambient COS concen- 
trations (Seibt etal. 2010), it is reasonable to assume that 
the COS mole fraction at the CA reaction site (j£ s ) is much 
smaller than its ambient concentration. We therefore follow 
others (Seibt etal. 2010; Stimler etal. 2010a) in assuming 



that can be neglected, as has been done on the right-hand 
side of Eqn 3. As depicted in Fig. 1, the similarity in diffu- 
sion pathways of C0 2 and COS depends on where actually 
most of COS becomes hydrated by CA, which is somewhere 
between the plasma membrane and the chloroplast stroma. 
Additional differences in gi and gf arise from differing 
biochemical reaction rates of CA and Rubisco, respectively, 
which are implicit in their numerical values. 

In summary, the diffusion pathway of COS is not identical 
to, but more similar for C0 2 than Ff 2 0 (Fig. 1). In the fol- 
lowing, we will discuss the implications of these differences 
for using COS as a tracer for canopy C0 2 and H 2 0 
exchange. As previous studies were almost exclusively inter- 
ested in linking photosynthesis to COS exchange, we first 
develop the procedure for estimating C0 2 from COS fluxes 
and then turn to the subject of water vapour fluxes. We 
would like to note that the approach outlined below, in 
contrast to what several studies have suggested (Stimler 
et al. 2010a), allows the quantification of canopy net, but not 
gross photosynthesis. Gross photosynthesis equals net pho- 
tosynthesis minus any autotrophic respiration that contin- 
ues in the light (Larcher 2001), a dissimilatory process with 
no apparent link to leaf COS uptake. 



ESTIMATING CANOPY NET PHOTOSYNTHESIS 
FROM COS EXCHANGE MEASUREMENTS 

Linking leaf-level COS and C0 2 exchange 

In order to provide independent estimates of net photosyn- 
thesis, leaf COS exchange measurements have to allow 
eliminating the unknowns in Eqn 1. Given that Eqn 3 con- 
tains three unknowns (g b s , gs S and gj S assuming fi s and % a s to 
be known) and inserting Eqn 3 into Eqn 1, making use of 
Eqn 4, would add two new unknowns for the one removed, 
this is however not possible. This problem has been over- 
come by parameterization of the relationship between COS 
and C0 2 fluxes through the ratio of their deposition veloci- 
ties, that is, the flux, normalized with the ambient concen- 
tration, of COS relative to C0 2 (for negative values of Ff), 
that is 



Xa 



(5) 



Based on an estimate of X c and measurements of ^ a c , 
and F] S , Ff (negative values only) may then be simply cal- 
culated as (Campbell et al. 2008): 



X^Fl_ 



(6) 



Up to now, parameterizations of X c have treated it as a 
vegetation-type specific constant (Sandoval-Soto et al. 2005; 
Campbell et al. 2008; Seibt et al. 2010), but how constant is 
X c expected to be? In order to answer this question, we 
rearrange Eqns 1 and 3 to yield the respective deposition 
velocities, that is 
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Making use of the diffusivity ratios (Eqn 4), we now intro- 
duce a non-dimensional coefficient j3° as 
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which when combined with Eqn 5 yields the following 
expression: 
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A similar, albeit less general, expression was derived by 
Seibt et al. (2010). It functionally relates the deposition 
velocities for CO2 and COS and shows that the CO2 depo- 
sition velocity will, for any given COS deposition velocity, 
vary with four unknowns: the ratio of intercellular to 
ambient CO2 concentration [cf. fig. 3 in Seibt et al. (2010)], 
gb S , g s s and gi S , which have been incorporated into our coef- 
ficient ff for convenience. 

By expressing these three conductances in a mutual 

„ s _ s 



fashion as ratios, that is, Rw, = 1 



and R,j = : 



, and after 



some manipulation of Eqn 9, the roles of g b s , g s s and gi S in 
determining pF can be explored, that is 
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Two limits of Eqn 11 are useful to explore: the first one 
relates to the case when R hs approaches infinity, as might be 
observed in a well-ventilated leaf enclosure where g b » g s . 
In this case, it can be shown that = 1.21 (1 + R si ), that is, pF 
increases in a linear fashion with the ratio of stomatal to 
mesophyll conductance with the slope and y-intercept equal 
to 1.21, that is, the g s s /g s c ratio. The second limit relates to 
the case of an infinite gf, that is, R si approaches zero. In this 

1 

^bs 



1 



case, Eqn 11 reduces to /3 



1 

1.21 



which corre- 



1.14fl bs 

sponds to an asymptotic increase from 1.14, that is, the 
gb S /gb C ratio, when _R bs is very small to the previously found 
limit of 1.21 when R bs becomes very large. Assuming 



physically/physiologically plausible values of 1-100 and 
0.001-2 for R bs and R si , respectively, yields pF values in the 
range of 1.2-3.6 (Fig. 2). Combining these values with %Fl 
X<F ratios typical for C3 plants of 0.5-0.8 (Larcher 2001) 
yields a range of 0.6-4.3 for A c , broadly in correspondence 
with the spread of 0.4-10.3 reported by Sandoval-Soto et al. 
(2005) in a recent literature survey. Lower %Fl%Jr ratios of 
C 4 plants result in correspondingly lower A c values (Stimler 
et al. 2011). An example of the actual variability of X c as a 
function of the intercellular to ambient CO2 concentration 
and pF is shown in Fig. 3 using data digitized from fig. 6 of 
Stimler et al. (2010a). While part of the variability in X° 
relates to changes in the intercellular to ambient CO2 con- 
centration ratio (leading to the asymptotic decrease as %f/ 
% a c decreases), it is also clear that part of the between- and 
within-species variability has to be ascribed to variability in 
pF and associated changes in R bs and R si . 

An example of how diurnal variability in g b s and g s s may 
affect pF, and thus in turn A c , under field conditions, is given 
in Fig. 4 using boundary layer and stomatal conductance 
values measured over the course of 1 d in a temperate 
mountain grassland (Wohlfahrt et al. 2010). In this particu- 
lar case, differences in pF values between the upper and 
lower canopies (1.4-1.8 and 1.2-1.3, respectively) were 
driven mainly by differences in stomatal conductance 
(assuming a constant g; s of 0.3 mol m~ 2 s -1 ), as Rbs values 
(6-21) were in a similar range in the upper and lower cano- 
pies. Assuming a %Fl%^ ratio of 0.7 (Larcher 2001) results in 
2F values of 1.8-2.3 and 2.5-2.6 in the upper and lower 
canopies, respectively. On the other hand, assuming a con- 
stant gi S is likely to be incorrect because experimental evi- 
dence suggests CA activity, which is implicit in gi S , to be 
dependent on leaf cellular pH (Stiltemeyer 1998), light 
intensity (Moskvin et al. 2000), as well as to be under circa- 
dian control (Eriksson et al. 1998). 

Another complicating issue, that up to now has not 
received appropriate attention, is whether A c determined 
from measurements in well-ventilated leaf enclosures is 
applicable to real-world, within-canopy transfer processes. 
Firstly, within plant canopies boundary layer conductances 
are not necessarily large as opposed to stomatal conduc- 
tances (Baldocchi 1988), for example, combination of low 
wind speeds and vigorous transpiration. Secondly, transfer 
mechanisms across the boundary layer may be dominated 
by free instead of forced convection with light winds 
and strong leaf-to-air temperature gradients (Finnigan & 
Raupach 1987), which would change the exponent in Eqn 4 
to 0.75 (Campbell & Norman 1998). 

Both the variability in 2F inferred from theoretical consid- 
erations (Fig. 2) and experimental evidence shown in Figs 3 
and 4 suggest A c to differ between plant species and vary 
dynamically in response to changing environmental condi- 
tions. This conclusion seems to be inconsistent with constant 
A c values in the range of 2-3 reported in previous studies 
(Sandoval-Soto et al. 2005; Seibt et al. 2010). On the basis of 
Eqns 10 and 11, we suggest 2F values on the order of 2-3 to 
reflect similar experimental conditions rather than an under- 
lying universal principle. Most of the studies to date have 
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Figure 2. Coefficient ff (Eqn 11) as a function of the relative magnitudes of boundary to stomatal (gb7g s s ) and stomatal to internal 
igflgf) conductance for COS. 



been performed at leaf level under laboratory conditions by 
using enclosure systems [see review by Sandoval-Soto et al. 
(2005)] and thus tend to be biased towards high boundary 
layer (due to leaf enclosures usually being well ventilated) 
and relatively low stomatal conductances (due to limitations 
resulting from pot size and growth under relatively low light 
conditions) . For example, for an infinite R hs and a xf/XsP ratio 
of 0.7, X c values between 2 and 3 are observed for R s i < 0.4. 
Limited support for this hypothesis derives from the recent 
work of Sandoval-Soto et al. (2005), where field as opposed 
to laboratory studies tended to yield higher X c values. If true, 
the reliability of available A c values for estimating canopy 
net photosynthesis under field conditions may be seriously 
questioned. 



Transferring the leaf-level COS-C0 2 
relationship to the canopy level 

The ecosystem level net exchange of COS (F e s ) comprises, 
in addition to the leaf exchange, flux contributions by 
above- and below-ground sources/sinks, that is 



Ff=F?+Fj+F*. (12) 

Here, FJ 1 represents above-ground COS flux contributions 
other than from leaves (i.e. woody organs, flowers, attached 
dead plant matter, etc.), and F s s represents COS fluxes 
from/to the soil surface. Both COS emission from and uptake 
by soils have been reported, although COS emissions from 
soils in some older studies appear to be due to experimental 
artefacts caused by the use of COS-free instead of ambient 
air in soil enclosures (Watts 2000). Recent studies, where soil 
chambers were flushed with ambient air, have reported soil 
COS uptake rates relative to deposition to vegetation that 
ranged from negligible (Xu et al. 2002; White et al. 2010) to 
dominant (Kuhn et al. 1999). Velmeke (1993), cited in Xu 
et al. (2002), investigated COS exchange in branches with 
and without leaves and found deposition and emission of 
COS, respectively. Given the scarcity of ecosystem-scale 
COS flux measurements (Xu et al. 2002), the significance of 
soil and non-leaf COS exchange is thus unclear and F n \ + F s s 
in Eqn 12 should not be neglected a priori. Ideally, F n i s + F s s 
would be quantified concurrently with F e s in order to derive 
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Figure 3. The ratio of COS to CO2 deposition velocities (A c ) as a function of the intercellular to ambient CO2 concentration 
(1 - Xi C /Xa C )- Lines represent simulations based on Eqns 9 and fO for various combinations of boundary to stomatal (Rbs) and stomatal to 
internal (R a i) COS conductance ratios (the resulting jf values are given in the right lower corner). Symbols represent data digitized from 
fig. 6 of Stimler et al. (2010a). Ambient CO2 and COS concentrations refer to 380 //mol moH and 500 pmol mol -1 , respectively. 
Calculations assumed a constant boundary layer conductance to CO2 (1.17 mol nT 2 s _1 ; fig. 3 in Stimler et al. 2010a). 



canopy-scale Fi by difference. Such an approach may 
however be problematic in practice, as there is a mismatch in 
footprint of soil/branch enclosure as opposed to above- 
canopy micrometeorological flux measurements, which may 
be aggravated by the presence of spatial heterogeneity 
in COS sources/sinks. In addition, concurrent measurements 
of 7v s + F s s and F e s would significantly increase the experi- 
mental effort. As an alternative, night-time measurements of 
F e s might be used for estimating daytime F^ s + F s s , similar to 
the current practice of estimating daytime ecosystem respi- 
ration from night-time C0 2 flux measurements (Wohlfahrt 
et al. 2005). For this approach to work, night-time leaf COS 
exchange should be negligible, as confirmed by Sandoval- 
Soto et al. (2005), who found COS uptake during darkness to 
virtually cease when stomata are closed. In contrast, White 
et al. (2010) reported significant COS uptake by loblolly pine 
during darkness. Micrometeorological night-time COS flux 
measurements may however be problematic due to method- 
ological limitations during calm and stable atmospheric con- 
ditions (Aubinet 2008). It remains to be determined whether 
approaches for dealing with unreliable night-time measure- 
ments developed for C0 2 , for example, filtering of data 



according to the magnitude of turbulence and imputation 
of resulting gaps based on empirical regression models 
(Goulden et al. 1996; Gu et al. 2005) are applicable to COS 
as well. 

As mentioned earlier, only a handful of concurrent 
ecosystem-scale COS and C0 2 flux measurements have 
been published so far (Hofmann, Hofmann & Kesselmeier 
1992; Bartell et al. 1993; Xu et al. 2002). In order to explore 
the magnitude and variability of ecosystem-scale X c and 
differences to leaf-scale values, we re-analysed COS, CO2 
and H 2 0 flux data collected by Xu etal. (2002) over a 
Norway Spruce forest in Germany using the relaxed 
eddy accumulation (REA) method. Due to the lack of 
concurrent measurements of F nl s + F s s and reliable 
night-time REA estimates of F e s , we have assumed F e s = Ff. 
This simplification is supported by soil COS exchange mea- 
sured at the same site by Steinbacher, Bingemer & Schmidt 
(2004) showing an average deposition of -0.81 ± 0.03 
pmol m" 2 s" 1 to the soil as opposed to an average ecosystem 
deposition of -93 ±11.7 pmol m~ 2 s _1 (Xu et al. 2002). Eco- 
system COS flux measurements of Xu et al. (2002), however, 
were on average positive during the morning and the 
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Figure 4. Diurnal time course of measured boundary layer (gb) and stomatal (g s ) conductance to COS and the resulting ff coefficient 
(assuming a constant gj s of 0.3 mol nr 2 s _1 ). Data are for Trifolium pmtense measured on 4 June 2009 in a temperate mountain grassland 
in Austria. Upper and lower canopies refer to 0.35 and 0.05 m canopy height, respectively. For further details regarding the experiment, 
we refer to Wohlfahrt et al. (2010). 



evening when the NEE was positive (i.e. net loss of CO2 to 
the atmosphere), suggesting that other COS sources (e.g. 
woody plant material) may be playing an important role. 
Disentangling canopy net photosynthesis from NEE mea- 
sured above the canopy would require estimates of soil and 
woody respiration, which were not available for this study. 
However, by following a common practice (Wohlfahrt et al. 
2005), we used night-time NEE measurements under windy 
conditions for parameterizing RECO as a function of air 
temperature. RECO was then extrapolated to daytime tem- 
perature conditions to derive GPP as GPP = NEE - RECO. 
We have used GPP as a proxy for canopy net photosynthesis, 
recognizing that canopy net photosynthesis is larger (i.e. less 
negative) than GPP by the amount of CO2 respired from 
autotrophic tissues during daylight (Larcher 2001). X c was 
then calculated separately for NEE and GPP by dividing 
the measured ecosystem-scale COS fluxes normalized with 
ambient COS mole fractions with the measured NEE 
(inferred GPP) normalized with ambient CO2 concentra- 
tions in analogy to Eqn 5. 

Values of A c calculated by considering NEE and GPP 
averaged 10.3 and 5.8 (Fig. 5), respectively. Considering that 
GPP overestimates canopy net photosynthesis, A c values 
resulted to be at least by a factor of 2 larger than the range 
of 2-3 reported by Sandoval-Soto etal. (2005) and Seibt 



et al. (2010) based on leaf-level laboratory enclosure 
studies. On a diurnal timescale, A c varied from 9 to 12.4 
(38%) and from 4.4 to 7.1 (59%) for NEE and GPP, respec- 
tively. The range of X c variability is comparable in magni- 
tude to the one deduced from diurnal changes in g b s and g s s 
(31%; Fig. 3), directly supporting our idea that changes in 
gb S , gs S (and most likely gi S ) cause diurnal variability in fF 
and thus in X c . With known A c values, Eqn 10 can be used 
to explore possible values for {f and X^X^- Commonly 
observed jfi C /^ a c ratios between 0.5 and 0.8 resulted in $f 
values of 0.35-0.87, well below the theoretical limit of 1.14 
calculated from Eqn 9. In order to reach ff values within 
the theoretical range (Fig. 2), xf/Xa C values of at least 0.85 
are required. While X'^X^ values > 0.8 are higher than the 
usual range (Larcher 2001), Huber (1993), cited in Seibt 
et al. (2010), indeed determined values of 0.91-0.93 

for Norway Spruce, resulting in of 1.9-2.5 (Fig. 5), that is, 
clearly within the range shown in Fig. 2. 

ESTIMATING CANOPY TRANSPIRATION FROM 
COS EXCHANGE MEASUREMENTS 

Although several authors have linked ecosystem H 2 0 and 
COS flux measurements (Hofmann etal. 1992; Xu etal. 
2002), to our knowledge only one study (Seibt et al. 2010) 
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Figure 5. Ecosystem-scale A c (a) and X v (b) values 
re-calculated from bin-averaged CO2, H2O and COS flux 
measurements over Norway Spruce (Xu et al. (2002); note that 
periods with average net COS emission were excluded). GPP was 
calculated by parameterizing RECO, derived from night-time 
NEE measurements under windy conditions (horizontal wind 
speed > 3 m s -1 ), as a function of air temperature which was then 
extrapolated to daytime conditions as GPP = NEE - RECO. jf 
values in (a) were calculated from measured X c values (based on 
GPP) and an # c /£ a c ratio of 0.92 [Huber (1993) cited in Seibt 
et al. (2010)]. Dashed lines in (b) refer to X v at various RHs 
calculated from ff shown in (a) and an average f^ljf ratio 
of 1.5. 



explicitly mentioned the potential of COS as a tracer for 
estimating canopy transpiration. This may be due to the fact 
that the diffusion pathway of COS has less in common with 
H 2 0 than CO2 (Fig. 1) and COS may thus be anticipated to 
be less suitable as a tracer for H 2 0. 

In analogy to Eqn 10, the following link can be estab- 
lished between transpiration and the COS deposition 
velocity: 



(13) 



The coefficient j3 v can be calculated in analogy to Eqn 9 but 
does not equal fp because of differences in the COS to 



C0 2 and H 2 0 diffusivity ratios. The ratio of [3 y /[f varies 
between ca. 1.4 and 1.7, depending on R bs . 

Provided leaf temperature and atmospheric pressure 
(required to estimate the saturation water vapour mole 
fraction %i V ) are available and assuming a known value for 
j3 v (our previous discussion on £F applies here as well), 
Eqn 13 can in principle be used to infer transpiration on the 
basis of a measured COS deposition velocity. While leaf 
temperature is usually known with enough precision in leaf 
enclosures [but see Sharkey, Wiberley & Donohue (2008)], 
this is not granted at canopy scale under field conditions. 
Therefore, we further develop Eqn 13 according to Penman 
(1948) by introducing 7\ = T a + AT, where T„ represents air 
temperature and AT the difference between T\ and T B . As a 
consequence, in Eqn 13 can now be reformulated as: 



Xi e { e s (7j) _ e s (TJ + sAT 
~ v „ ~ „ ' 



(14) 



which holds for small values of AT and where e represents 
the vapour pressure and s the slope of the saturation vapour 
pressure function (Campbell & Norman 1998). Combining 
Eqns 13 and 14 yields, after some re-arrangement, the fol- 
lowing expression: 



VPD + sAT 



(15) 



where VPD refers to the vapour pressure deficit of air, that 
is, e s (T a >e a . 

Equation 15 clearly addresses the relationship between 
transpiration and COS deposition velocity showing how 
that relationship may be sensibly affected by changes in AT, 
T„ and e a , independently of /f . If AT = 0 (isothermal case), 
the right-hand side of Eqn 15 reduces to ((RFr 1 - 1) /3 V )" 1 , 
scaling positively with RH. When AT # 0, additional non- 
linear variability proportional to AT arises because 5 
increases exponentially with temperature (Campbell & 
Norman 1998). 

In contrast to the ratio, which appears relatively 

conservative under a wide range of environmental condi- 
tions (Larcher 2001), the ratio of the saturation vapour 
pressure (at leaf temperature) to ambient vapour pressure 
may be expected to be much more variable, causing larger 
variability in X v as compared with X c . Radiometric leaf 
temperature measurements would greatly reduce this 
problem in principle, may though be difficult in practice due 
to contributions by non-transpiring components such as the 
soil surface and/or woody plant tissue (Sanchez et al. 2009). 
In addition, radiometric measurements of leaf temperature 
may not necessarily represent a good estimate to be used 
for calculating the intercellular saturation vapour pressure 
due to within-canopy differences in the sources of thermal 
radiation and latent heat (Campbell & Norman 1998). 

Because to our knowledge no data regarding X y have 
been published so far, we again make use of the study of Xu 
et al. (2002) to investigate the magnitude and diurnal vari- 
ability of X Y at ecosystem scale. As shown in Fig. 5, A v results 
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to have an average value of 1.64 (ranging from 1.2 to 2.2), 
exhibiting higher and lower values in the morning/evening 
and noon, respectively, as expected from the positive rela- 
tionship with RH and the typically lower RH values around 
noon. Assuming a constant ft* Iff ratio, the additional vari- 
ability in A v (81%) as opposed to GPP-based A c (59%) is 
due to diurnal variability in the X?l%* ratio, supporting our 
arguments raised previously. 

ESTIMATING CANOPY STOMATAL 
CONDUCTANCE FROM COS 
EXCHANGE MEASUREMENTS 

By re-arranging Eqn 3, leaf stomatal conductance to COS, 
which may be converted to CO2 or H?0 via Eqn 4, may be 
calculated as: 



1 



1 



(16) 



Equation (16) differs from the equations developed above 
for canopy net photosynthesis and transpiration in three 
different aspects: (1) the COS to C0 2 (A c ) and H 2 0 (A v ) 
deposition velocity ratios are not included; (2) instead, 
Eqn 16 contains two unknowns as absolute numbers - the 
Gb S and the Gi S ; (3) Gb S and Gi S represent bulk quantities: 
Gb S includes both the aerodynamic and the quasi-laminar 
boundary layer conductance (Monteith & Unsworth 1990), 
while Gi S is the integral of the internal conductance over the 
canopy leaf area - in order to emphasize this difference to 
the leaf-scale conductances used in Eqns 1-4, 7-11, we have 
used uppercase letters in Eqn 16. While models for the com- 
bined aerodynamic and quasi-laminar boundary layer con- 
ductance are available [but see Liu et al. (2007) for 
systematic uncertainties], little a priori knowledge is usually 
available on the magnitude of gi S , which is inter alia why A c 
and A v were introduced earlier, prohibiting a reliable 
up-scaling to Gi S . A better quantitative understanding of gi s 
(cf. Stimler et al. 2011) needs to be achieved to employ COS 
as a practical tool for estimating canopy stomatal conduc- 
tance. 



CONCLUSIONS 

The theoretical basis for the observed co-variation of leaf 
net photosynthesis, transpiration and COS uptake was 
reviewed. In addition, the assumptions which are required 
to use COS as a tracer for canopy net photosynthesis, tran- 
spiration and stomatal conductance, were discussed. Based 
on our analysis, we identified the following two priorities for 
future research: 

1 We showed that A c and A v are not constants, but vary as 
a function of xf/% a c and ;fc v /£a V , respectively, and the 
ratios of boundary layer to stomatal and stomatal to 
internal conductance (Eqns 9-10 & 13). Due to the com- 
parably more conservative nature ofxf/% B c and the more 
similar diffusion pathway, our theoretical analysis 



suggests COS to represent a better tracer for CO2 than 
H 2 0. For routine application of Eqn 6 in field conditions, 
it will be necessary to develop a better understanding of 
the in situ variability and the factors controlling A c and 
A v , which may result to be quite different from what has 
been observed under laboratory conditions by using leaf 
enclosures. To this end, it will be more important to study 
the four component processes controlling A c and A v in 
greater detail than limiting the investigation to the depo- 
sition velocities ratios. A particularly important step into 
this direction, which is prerequisite to improving the 
potential of COS as a tracer for canopy conductance, is to 
better understand and quantify variability in gj S between 
species and on diurnal and seasonal time scales. A prom- 
ising approach for independently characterizing A c under 
field conditions has already been put forward by Seibt 
etal. (2010), who showed that X c is related to the 13 C 
discrimination during photosynthesis and that for 
example 8 13 C of leaf samples could be used as time- 
integrated estimates for X c . 
2 The potential of COS as a tracer for canopy-scale 
exchange processes hinges upon our ability to operation- 
ally disentangle leaf from any other ecosystem COS 
exchange. While progress has been made in understand- 
ing the factors which drive soil COS exchange (Kuhn 
et al. 1999; Van Diest & Kesselmeier 2008), we definitely 
need a better quantitative understanding of the contribu- 
tion of soil and other non-leaf ecosystem components 
to the overall ecosystem-scale COS exchange. Here it 
appears that advantage can be taken from the vast 
amount of experience gathered in recent years in disen- 
tangling ecosystem CO2 flux components (Reichstein 
et a/. 2005). 

In summary, our study confirmed previous pioneering work 
highlighting the potential of COS as a tracer for canopy net 
photosynthesis; however, we have also indicated a number 
of limitations. For the first time, we have assessed the link 
between leaf COS and H 2 0 exchange, addressing the 
potential and the limitations of COS as a tracer for canopy 
transpiration and stomatal conductance. At present, we face 
a serious lack of ecosystem-scale field measurements (Briihl 
et al. 2011) that represent an essential requirement for 
assessing the practical significance of these limitations and 
whether or not ecosystem COS flux measurements will be 
able to provide sensible constraints on canopy net CO2 and 
H 2 0 exchanges. Such measurements will be of great value 
also for studies at larger scales (e.g. regional or global), 
which aim at inverting concurrent COS and CO2 concen- 
tration measurements to disentangle the components of the 
net ecosystem CO2 exchange (Campbell etal. 2008). We 
thus advocate concurrent measurements of ecosystem-scale 
COS, CO2 and H 2 0 exchange and the corresponding com- 
ponent fluxes to allow testing and validation of our theo- 
retical COS exchange models and their relation to canopy 
C0 2 and H 2 0 fluxes under field conditions. Here, great 
advancements can be expected from recent developments 
in analytical instrumentation (Graus etal. 2010; Stimler 
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et al. 2010b) that will allow quantifying ecosystem-scale 
COS exchange by using the eddy covariance method in the 
near future. 
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